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^ i' Abstract. We present results from numerical Classical Yang Mills calculations of 

' the initial stage Glasma field configurations in a relativistic heavy ion collision. 

We compute the initial gluon multiplicity from RHIC to LHC energies. The 
initial conditions for the classical field computation are taken from a dipole model 
parametrization tested on HERA data, meaning that all the parameters, including the 
normalization and the value of the saturation scale, are fixed from DIS data. 
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1> . PACS numbers: 24.85. +p, 25.75.-q, 12.38.Mh 

1. Introduction: the Glasma 

> 
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^ I At high energy the small x part of the wavefunction of a hadron or nucleus can be 

^ I understood as a classical color field radiated by static color sources formed by the large 

I X degrees of freedom [1]. This description, known as the color glass condensate (for 

O ' reviews see e.g. [21 [3]), provides a common framework for understanding both small x 
00 ' 
O 



X 



deep inelastic scattering (DIS) and the initial stages of relativistic heavy ion collisions. 
The cross section for small x DIS can be expressed in terms of the correlator of two 
Wilson lines (i.e. the dipole cross section), and the initial condition for the classical 



^ ' fields that dominate the first fraction of a fermi of a heavy ion collision is determined 
by these same Wilson lines. 

The earliest stage of an ultrarelativistic heavy ion collision is a coherent, classical 
field configuration of two colliding sheets of Color Glass. This first fraction of a fermi of 
the collision, in transition towards a equilibrated quark gluon plasma, is what we refer 
to as the Glasma [H [5] . The color fields of the two nuclei are transverse electric and 
magnetic fields on the light cone. The glasma fields left over in the region between the 
two nuclei after the collision at times < r < l/Qs are, however, longitudinal along the 
beam axis. They fields depend on the transverse coordinate on a length scale of order 
1/Qs- As the system expands the fields are diluted and can be treated as particles. This 
lowest order contribution, in the MV model, has been computed numerically by Classical 
Yang-Mills (CYM) computations of gluon production in heavy ion collisions [HI [71 [8]. 

In order to relate calculations of the initial gluon multiplicity to final state 
observables at RHIC and LHC, some simplifying assumptions have to be made. We will 
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Figure 1. Left : The Wilson line correlator vs. /g^lJ- in the MV model with different 
diseretizations of the longitudinal coordinate. The value Ny = 1 corresponds to the 
implementation in the CYM calculations. Right: The saturation scale (in the adjoint 
representation) in units of ^ for different values of the infrared regulator m used in 
the calculation. The straight and dashed lines correspond to the analytical {Ny oo) 



estimate 
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with different values for the constant C. The value 



corresponding to the CYM calculations is Qs = 0.57p /i (at m = and Ny ~ 1) 



in this simple study assume that the initial gluon multiplicity is approximately equal to 
that of the final state hadrons (charged and neutral); about 1000/1150 per unit rapidity 
at RHIC with = 130^4 GeV/200A GeV. This assumption can be justified by either 
rapid thermalization of the system and entropy conservation during its hydrodynamical 
evolution, or the phenomenological assumption of parton-hadron duality between the 
gluon and hadron multiplicities. Note that for the relation between the initial and final 
energy densities these two assumptions lead to very different interpretations. 

2. Value of the saturation scale 

To make genuine predictions for the multiplicities in heavy ion collisions we must be able 
to obtain the value of the saturation scale from DIS measurements at HERA. For the 
purposes of this study we will do this using two dipole cross section parametrisations, 
where the impact parameter dependence is modeled and constrained using diffractive 
HERA data. These are the IPsat (or Kowalski-Teaney) and the bCGC parametrisations. 
The IPSat model [9l|T0], is built from an eikonalised DGLAP-evolved gluon distribution 
function. Thus it contains the right color transparency limit at large . The bCGC 
model [H] is an impact parameter dependent version of the IIM [12] parametrization, 
which includes the main effects of BK evolution. For RHIC energies the values of the 
saturation scale at median h from these parametrizations are are Qs ~ 1-1 GeV (IPsat) 
and Qs ~ 1.0 GeV (bCGC) [lO]. 

So far the CYM computations of the Glasma initial state have been performed with 
initial conditions from the MV model [1], where the color charge density is determined 
by a phenomenological parameter g^[i. In order to use DIS data to set the value for 

one must consistently determine the value of the Wilson line correlation length 
(the saturation scale) in implementation of the MV model used in the numerical 



The glasma initial state at the LHC 



3 



o 




Figure 2. Left: Momentum space Wilson line correlator (scaled by k^"') in the IPsat 
model at RHIC energies ( "goal" ) and the same correlator from two field configurations 
constructed to reproduce it, on average. Center: Spectrum of produced gluons at LHC 
energy using IPsat and bCGC dipole cross sections. Right: Energy dependence of 
charged particle multiplicity from the IPsat and bCGC models, from the MV model 
assuming ~ and from a A + i? Ins-fit to lower energy data. 
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Table 1. Simrmary of CYM results. Calculated on a 5122-latticc at r = l.G/Qs- 

computations. This comparison has only recently been performed in Ref. [13] (see 
also the related discussion in Ref. [H]). The results are summarized in Fig. [H The 
obtained value Q^jg^ii ~ 0.57 means that the HERA-extracted values quoted above 
are in remarkable agreement with the earlier estimate g^[i ^ 2 GeV based on RHIC 
phenomenology alone [H [8] . 

3. Beyond the MV model, LHC predictions 

It is possible to go beyond the MV model in the CYM calculation of the glasma fields. 
This can be done by supplementing the given Wilson line correlator with additional 
assumptions of the statistics of the large x-sources. We adopt the commonly used 
"nonlocal Gaussian" approximation, in which the color fields y4"''(xT,x~) in the Wilson 
line U = P exp{ig J dx'A'^} at different values of the longitudinal coordinate x~ 
are taken to be independent and have a Gaussian distribution. With this additional 
assumption one can numerically construct the Wilson lines; the result of this procedure 
is shown in Fig. [2] (left). At large momenta the k^-space correlator is negative (absolute 
value plotted); only the positive part is used. This leads to a problematic behavior of 
the gluon spectrum at large (see Fig. [Reenter), but for bulk quantities, such as the 
total multiplicity, the effect is not significant. 

It is convenient to parametrize our results in terms of the "gluon liberation 
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coefficient" c [T3] which is defined using the ratio of the gluon muhiphcity to by 
,■2^^, = c^^^. The MV model calculation leads |13| to the result c ~ 1.1. We 
find that both the IPsat and bCGC models lead to very similar values, within 10% 
of the MV model result, for the liberation coefficient c. The gluon spectrum obtained 
using the IPsat or bCGC parametrisations is shown in Fig. [2] (center). There is a 
remaining uncertainty in the value of x used to determine the saturation scale in an 
ion-ion collision where, unlike in DIS, it is not precisely determined by the collision 
kinematics. Parametrically one should take x = c^Qs/ \/s with of order one. Taking 
= 1 yields a lightly smaller gluon multiplicity than our assumption based on the 
measured RHIC charged multiplicity. In making a prediction for LHC we adjust the 
constant in order to obtain dNg/dt] = 1150 at RHIC 200 GeV. Our first preliminary 
results for RHIC and LHC energies from this study are summarized in Table [1] 

The a/s dependence of the charged particle multiplicity from RHIC to LHC 
energies is shown in Fig. [2] (right). There is a significant difference between the 
bCGC and IPsat models, resulting from the much weaker energy dependence of Qs 
in the former. Thus the LHC measurements of particle multiplicities will significantly 
constrain our understanding of the high energy evolution in the Color Glass Condensate. 
In the meanwhile, to make LHC prediction more accurate, further theoretical and 
phenomenological work is needed to understand effects such as the impact parameter 
dependence of the dipole cross section and running coupling and pomeron loop effects 
on the evolution. 
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